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Q\ • Using the mode-by-mode summation technique the zero point energy of the 

, electromagnetic field is calculated for the boundary conditions given on the 

i , surface of an infinite solid cylinder. It is assumed that the dielectric and 

§ H | magnetic characteristics of the material which makes up the cylinder (e±, fi±) 

i-^h ■ and of that which makes up the surroundings (e2, ^2) obey the relation £\^i = 

■ £2/^2- With this assumption all the divergences cancel. The divergences are 

regulated by making use of zeta function techniques. Numerical calculations 
are carried out for a dilute dielectric-diamagnetic cylinder and for a perfectly 
conducting cylindrical shell. The Casimir energy in the first case vanishes, 
and in the second is in a complete agreement with that obtained by DeRaad 
and Milton who employed a Green's function technique with an ultraviolet 
regulator. 
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I. INTRODUCTION 



Calculation of the Casimir energy for nontrivial boundary conditions is a subject of 
intense ongoing activity In spite of this, the experience accumulated in this area still does 
not allow one to predict, without involved calculation, even the sign of this energy In 
this connection the development of new effective methods of calculating the Casimir energy 
is doubtless of interest. 

In recent papers HQ the zero point energy of the electromagnetic field with boundary 
conditions given on a sphere has been calculated by a method of direct summation of the 
eigenfrequencies. A substantial point here was the use of contour integration in the complex 
frequency plane ||. The divergences in this problem were removed by two subsequent 
steps. The first was accomplished by the subtraction of the vacuum energy of an infinite 
homogeneous space and the second made use of the technique of zeta function regularization. 
The same method was also applied for calculating the Casimir energy of a massless scalar field 
obeying Dirichlet and Neumann boundary conditions on a sphere ||. The same technique 
has recently been used to rederive the electromagnetic Casimir energy for a spherical shell || . 
The possibility of incorporating in this approach the dielectric and magnetic properties of the 
media was demonstrated in Q. These calculations are distinguished by a certain concision 
and simplicity in comparison with the Green's function techniques 0, which, however, retain 
superiority with respect to physical interpretation. 

The present paper seeks to show the efficiency of the mode summation method in calcu- 
lating the zero point energy of an electromagnetic field when the boundary conditions are 
given on an infinite circular cylinder. When applying the Green's function techniques [|],[7||§, 
this problem turns out to be more complicated than the corresponding calculation of the 
Casimir energy for sphere [p| P~5f| . 

The layout of the paper is as follows. In Sec. II the general integral representation 
is derived for the Casimir energy of an infinite solid cylinder surrounded by an uniform 
medium. The permittivity and permeability of the cylinder material (e\, /xi) and those of 
the surroundings (£2, ^2) are considered to be arbitrary. In principle they may depend on 
the frequency of the electromagnetic oscillations (dispersion of the media), but this point is 
beyond the scope of the present paper. In Sec. Ill it is assumed that the electromagnetic 
characteristics obey the condition eifii = 62^2 = c~ 2 where c is the speed of the light 
in the media (in units of that in the vacuum). When this condition is satisfied, all the 
divergences cancel between interior and exterior modes. Those divergences are regulated by 
employing the zeta function technique. In Sec. IV the cases when ^ 2 C 1 and £ 2 = 1 are 
considered numerically, £ 2 being (e\ — e 2 ) 2 / (£i + S2) 2 - The first case gives the Casimir energy 
of a dilute dielectric-diamagnetic cylinder, while the second case corresponds to a perfectly 
conducting infinitely thin cylindrical shell. Remarkably, the Casimir energy obtained for a 
tenuous medium vanishes, as it does for a tenuous dielectric cylinder. The result obtained in 
the second case is identical to that obtained by the Green's function method of calculating 
the energy and regulating the divergences by use of an ultraviolet regulator M. In Sec. V 
(Conclusion) the significance of the universal results obtained are discussed. 
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II. INTEGRAL REPRESENTATION FOR THE CASIMIR ENERGY 



We shall consider the following configuration. An infinite cylinder of radius a is placed 
in an uniform unbounded medium. The permittivity and the permeability of the material 
making up the cylinder are E\ and /^i, respectively, and those for surrounding medium are 
e 2 and /x 2 - It is assumed that the conductivity in both the media is zero. We will compute 
the Casimir energy per unit length of the cylinder. 

In the mode summation method, the renormalized Casimir energy is defined by 

z to 

where uj p are the classical eigenfrequencies of the electromagnetic oscillations in the system 
under consideration, and uj p are those in the absence of any boundary, that is, when either 
medium fills all space. (When the precise meaning is not required, we denote this by the 
formal limit a — > oo.) The set {p} stands for a complete set of quantum numbers (discrete 
and continuous) which is determined by the symmetry of the problem. Either sum in Eq. 
( P-1P diverges, therefore a preliminary regularization is required. 

In order for the eigenfrequencies to be found one needs to solve Maxwell's equations 
for the given configuration with allowance for the appropriate boundary conditions on the 
lateral surface of the cylinder. As is well known, it is sufficient to require the continuity 
of the tangential components of the electric field E and of the magnetic field H ||16|| . In 
terms of the cylindrical coordinates (r, 9, z) the eigenfunctions of the given boundary value 
problem contain the multiplier 

exp(—iu)t + ik z z + in6), (2-2) 

and their dependence on r is described by the cylindrical Bessel functions J n for r < a and 
by the Hankel functions of the first kind H n = for r > a. The eigenfrequencies are the 
roots of the equation (Ref. |TB[, p. 526) 

f n (k z ,LO,a) = 0, (2.3) 

where 

f n = A^A^^a, X 2 a) A™(A ia , X 2 a) - nWkl^fn - e 2f i 2 ) 2 ( J n (A ia ) H n {X 2 a)f , (2.4) 
with 

A™(Aia, X 2 a) = a/iiA 2 J^(Aia) H n (X 2 a) - a/i 2 Ai J n (Xia) H' n (X 2 a), 

A™(Aia, A 2 a) = aeiX 2 ^(A x a) H n (X 2 a) - ae 2 Xx J n (Xia) H' n (X 2 a), (2.5) 



x l = k i- k l k 2 i =e i ^ i u' 2 , i = l,2, n = 0,±1,±2,.... 

The indices TE and TM will be explained below. The prime on the functions J n and H n 
means differentiation with respect to their arguments. For given k z and n Eq. (|2.3|) has an 
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infinite sequence of roots u nm (k z ), m = 1, 2, . . ., these frequencies being the same inside and 



outside the cylinder [17]]. In view of this the Casimir energy (|2.1| ) can be rewritten as 

1 r°° dK 00 00 

n=— oo m=l 



1 f°° dk z 00 00 

£=-/ — 1 51 51 km(fc 2 ) - Wnm(fe)] , (2-6) 
Z J —oo ZTT __ , , __i 



where u) nm (k z ) stands for the uniform medium subtraction referred to above. 

The next step in our consideration is a representation of the sum in Eq. ( j2.6p in terms 
of the contour integral pfl 



1 /-oo &k ™ 11/" , , f n (k z ,uj,a) 
E = - —Y i u d u In '7, ' ' \ . 2.7 

Integration in ( |2.7| ) is conducted along a closed path C in the complex u plane which con- 
sists of two parts: C + which encloses the positive roots of Eq. ( |2.3|) in a counterclockwise 
sense, and C_ which encircles the negative roots in a clockwise sense. Therefore, we face 
to task of investigating the properties of the function f n (k z ,u,a) specifying the frequency 
equation. (However, the result ( |2.7| ) may be easily shown to be equivalent to the corre- 
sponding Green's function formulation. See Appendix A.) Generally this is a problem of 
extreme difficulty Therefore, in the next sections we shall consider specific cases introducing 
simplifying assumptions. 

The method of calculation of the Casimir energy proposed above can straightforward 
be generalized to the dispersive media. To this end, it is sufficient to treat the parameters 
Si and fii, i = 1,2 in the frequency equation (|2.3|) as given functions of the frequency uj. 
However, this issue is beyond the scope of the present paper. 



III. CASIMIR ENERGY OF AN INFINITE CYLINDER WHEN e x \t\ = £2/^2 

We assume that the permittivity and permeability of the cylinder material (e i and of 
the surroundings (E2, ^2) ar e not arbitrary but satisfy the condition 

siHi = e 2 fi2 = c" 2 , (3.1) 

where c is the light speed in either medium (in units of the speed of light in vacuum). The 
physical implications of this condition can be found in [T8^p0f . When equation ( |3 . 1| ) holds, 
we have Ai = A2 = A, and the frequency equation ( ^73|) is simplified considerably. It breaks 
down into two equations: for the transverse-electric (TE) oscillations 

A™(A, a) = Xa \pL X J' n {\a)H n (Xa) - n 2 J n (\a)H' n (\a)} = (3.2) 

and for the transverse-magnetic (TM) oscillations 

A™(A, a) = Xa [e x J' n {\a)H n {\a) - e 2 J n (Xa)H' n (Xa)} = 0. (3.3) 

In the general case [see Eq. ( |2.4| )1 such a decomposition occurs only for oscillations with 
n = 0. In Eqs. (|3.2|) and ( |3.3|) A is the eigenvalue of the corresponding transverse [membrane- 
like] boundary value problem |21[ 
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\ 2 = ~--k 2 z . (3.4) 



c 



Classification of the solutions of Maxwell's equations without sources in terms of the 



TE- and TM- modes originates in waveguide theory [|Tj],^T|,^] . The main distinction of the 
propagation of electromagnetic waves in waveguides, in contrast to the same process in 
unbounded space, is that a purely transverse wave cannot propagate in a waveguide. The 
wave in waveguide must necessarily contain either longitudinal electric or magnetic fields. 
The first case is referred to as the waves of electric type [transverse-magnetic (TM) waves] 
and in the second case one is dealing with waves of magnetic type [or transverse-electric (TE) 
waves]. This classification proves to be convenient in studies of electromagnetic oscillations 
in closed resonators as well. 

Replacing the function f n (k z ,u,a) in Eq. (|2.7|) by the left hand sides of Eq. ( |3.2| ) and 
( p.3|) and changing to the integration variable A we arrive at the following representation for 
the Casimir energy 

c r dk ~ 1 1 r fW-77,. . Al E (\a)A™(Xa) 

Here we have distorted the contour of integration to C = C' + + C'_. We take C' + to consist 
of a straight line parallel to, and just to the right of, the imaginary axis (—zoo, +ioc) closed 
by a semicircle of an infinitely large radius in the right half-plane. C'_ similarly is a line 
parallel to, and just to the left of, the imaginary axis, closed by an infinite semicircle in the 
left hand plane. On both semicircles the argument of the logarithm function in Eq. ( p.5| ) 
tends to 1. As a result these parts of the contour C do not give any contribution to the 
Casimir energy E. When integrating along the imaginary axis we chose the branch line of 
the function <f(X) = \J\ 2 + k 2 to run between — ik z and ik z , where k z = + \fk 2 > 0. In 
terms of y = Im A we have 





1.2 


A, 2 


-y 2 , 



y>k z 



tpjiy) = { ± s j kl-y l , \y\ < k z , (3.6) 
-ijy 2 - k 2 z , y < -k z , 



where the sign on the middle form depends on whether we are to the right or the left of 
the cut. Thus contributions to Eq. ( |3.5|) due to the integration along the segment of the 
imaginary axis (—ik z , ik z ) cancel between C' + and C'_, and Eq. (|3.5D acquires the form 



n=— oo 



A™ ( , 00 )atm ( , oo) • 



Changing the order of integration of k z and y and taking into account the value of the 
integral 

V dk z sjy 2 -kl = jy\ (3.8) 
o v 4 

we obtain after substitution ay — >• y 
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E 



E 



d v In • 



A™(iy)A™(n,) 



" '"ATE(, 00 ) A TM (z00; 

Further we shall need the modified Bessel functions I n {y) and K n (y) p3| 

J n (y) = i~ n J n {iy) 



K n {y) 



7T 



H n (iy) 



and their asymptotics at fixed n and y — > oo 



/n(j/) 



(3.9) 



(3.10) 
(3.11) 



(3.12) 



7T 



<(v) 



7T 



2y ' — V 2 ?/ 

With the help of this we derive from Eqs. (|3.2|) and (|3.3|) 
A™(%) 2y 



A™(ioo) /ii+/i 2 



[V\l'n{y) K n{y) - H 2 In(y)K' n (y)}, 

e 1 I' n (y)K n (y)-e 2 I n (y)K' n - 



A™(ioo) e 1 + e 2 
Making use of all this, we can recast Eq. ( |3.9|) into the form 



E 



Ana 2 



E 



ii U + e 1 + 2 



(e + e~ l ) Uy)I' n {y)K n {y)K' n {y) 



(3.13) 



(3.14) 



(3.15) 



Here a new notation e = S\/e 2 has been introduced, /x has been eliminated by condition 
( p.l|) , and when going from ( |3.9|) to fl3.15|) an integration by parts has been done, the 
boundary terms being omitted. The last point will be justified further when the removal 
of the divergences is discussed. The argument of the logarithm in Eq. (|3.15 ) is simplified 
considerably if one uses the value of the Wronskian of the modified Bessel functions I n {y) 
and K n (y) @ 



I n (y)K' n (y) - r n (y)K n (y) 



1 

y 



and the identity 

I n (y)K' n (y) + I' n {y)K n {y) = (I n (y)K n (y))' . 
Finally Eq. ( [3.15D acquires the form 



(3.16) 



(3.17) 
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E= £ E n , (3.18) 

n=— oo 

where 

En= 4^l y d y^{ l -?Wn{y)K n {y))'} 2 } (3.19) 
with £ = (1 — £)/(! This is a simple unregulated generalization of (4.5) of Ref. 0, and 



the cylindrical analog of the spherical form (7.1) of Ref. [12 



From the asymptotics (|3.12j) and ( [3.13] ) it follows that the integral in ( |3.19 ) diverges 
logarithmically when y — ► 00. At the same time the sum over n in Eq. (|3.18 ) also di- 
verges because at large n the uniform asymptotic expansion of the modified Bessel functions 
gives ]23| (see (4.6) of M) 



E, 



n— >oo 



- 1 Z =E°°. (3.20) 



16na 2 Jo (1 + z 2 ) 3 



Here the change of variables y = nz has been done. Disregarding for the moment that the in- 
tegral in Eq. (|3.20j ) is divergent, we employ here the Riemann zeta function technique |]24]j25 
for attributing a finite value to the sum in ( |3.18|) 

00 

E= (En-E^ + E 00 ) 

n=—oo 

00 00 

= £ {E n - E°°) + £ E™ 

n=— 00 n=— 00 

00 00 

= £ E n + E°° Y, n °> ( 3 - 21 ) 

n=— 00 n=— 00 

where _E n stands for the "renormalized" partial Casimir energy 

E n = E n — E°°, n = 0,±l,.... (3.22) 

We now have to treat the product of two divergent expressions E°° •£)J£ = _ 00 n ° m ore precisely, 
by presenting it in the following form 

E°° ■ Y n° = ^— lim / , Z a * • [2C(s) + 1] 

167ra 2 7o (1 + z 2 ) 3 L SV ! 1 

= —^2 H ^ f 1 - tV PCW «] (3-23) 



167ra 2 s^o+ Vs 4 

ce ,. /1 3 



167ra 2 s"-^o+ Vs 4 
c£ 2 



lim ■ [-ln(2vr)s] 



■ln(27r). 



167ra 2 

Finally, the Casimir energy acquires the form 
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E 



,_ _ lbna z 



(3.24) 



Now we deduce from Eqs. flT22p , flUSQ , and ( ggg ) for E, 



Ana 2 



ydyln 1 -f crjj/) 



+ +- 



£ 2 Z" 00 z 5 Gb 



4 7o l + 



2^3 



(3.25) 



n 



0,1,2, 



(3.26) 



where cr n (y) = y(I n (y)K n (y))' . Since both integrals in Eq. ( |3.25D rdiverge, the finite sum is 
to interpreted in a precise manner as specified below. 

The removal of the divergences by making use of the ( function justifies dropping the 
boundary terms in the integration by parts when we went from Eq. 



to (FT51) . To 

see this one can at first remove the divergences in Eq. (|3.9| ) by employing the Riemann 
C function as described above. After that the integration by parts in the formula for the 
renormalized partial Casimir energy E n can be done, with boundary terms now vanishing 
rigorously. Proceeding in that way one again arrives at Eqs. ( |3.25| ) and ( |3.26| ). 

In order for Eq. ( |3.25 ) to be cast in a form suitable for numerical evaluations both terms 
there should be placed under s single integral sign. To this end, for n / the change of 
variables z = y/n has to be done in the second term. This change is inverse to the one 
which has led to the asymptotic form (|3.20|) . This yields 



E n 



Ana 2 



+ 



e 



A (n 2 + y 



2\3 



For n = Eq. (|3.25|) is rewritten as 



E 



Ana 2 Jo 



ydy In 1 - £V 2 (y) 



e 



n 



1,2, 



A 1 + y 



2\3 



(3.27) 



(3.28) 



The integrals in these formula? converge because for y — > and we have [23 



In this limit 



I n {y)K n {y)^—. 



(3.29) 



(3.30) 



On the other hand, for large y 



(3.31) 



By making use of uniform asymptotics of the Bessel functions |23f we deduce from (|3.27|) 
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E, 



■oo 4ira 2 
Thus the Casimir energy 



10 - 3£ 2 28224 - 7344£ 2 + 720£ 4 Q ( 1 



960n 2 



15482880n 4 



E = E + 2 E n + 



ce 



n=l 



ln(27r) 



(3.32) 



(3.33) 



with E and E n , n > 1, given in ( |3.28| ) and ( |3.27| ), respectively, is finite f2~6| . One can 
advance further only by considering special cases and applying numerical calculations. 



IV. DILUTE COMPACT CYLINDER AND PERFECTLY CONDUCTING 

CYLINDRICAL SHELL 

We begin by addressing the case when £ 2 <C 1. Because we are assuming the condition 
E\\i\ = e 2 fii = c~ 2 and ^ 2 C 1, this is not the same situation as a dilute compact cylinder 
with |ei — e 2 | <C 1 and ^1 — ^2 — 1- In this case 

2 (si - e 2 f 



£1 ~ £2 



(4.1) 



where e — (e\ + E^fl. Retaining in Eq. ( |3.28|) only the terms proportional to £ 2 we obtain 



En 



ce_ 

Area 2 Jo 
4vra 2 1 



4(1 + y 2 ) 3 
0.490878) . 



To estimate £J n , n > 0, we can use the leading asymptotic behavior ( |3.32| ) 

oi 2 ( 1 7 



4vra 2 V96n 2 3840n 4 



(4.2) 



(4.3) 



To a precision of 10~ 6 , we evaluate Eq. ( |3.33| ) by substituting in the value of E , ( |4. 2| ) 
integrating E n numerically for n = 1, ... 5, and asymptotically using Eq. ( |3.32| ) for n > 6 
with the resultf] 

ce 2 ( 5 

s - I -0.490878 • :>\ 1 . 

48 



E 



Ana 2 

ce 

Ana 2 

ce 

Ana 2 



— , 1 00 1 7 00 1 1 
-0.490878 + + V -r + - ln(2 

^ n /IS ^ r,2 1Q9D ^„4 4 v 



n=1 ™ n=6 n- 1920 — 



7T 



(-0.490878 + 0.027638 + 0.003778 - 0.000007 + 0.459469) 
(0.000000) . 



(4.4) 



x The cancellations here are very severe. If the asymptotic approximation were used for all n, a 
positive result would be found, E ~ (c£ 2 /4-7ra 2 )(— 0.00108). Unlike for the spherical case, doing 
the integral exactly is essential. 
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Thus the Casimir energy of a cylinder possessing the same speed of light inside and outside 
proves to be zero! This is to be contrasted with the positive Casimir energy found for a 
dilute ball with the same property |I^ , |T2| ,P| , 

£baii « tJ-£ 2 = 0.046875^. (4.5) 
64a a 

It is further remarkable that the same zero result is found for a dilute dielectric cylinder, 
that is, one with \i = 1 everywhere and e > 1 inside the cylinder, and e — 1 outside, a 
result which may be most easily confirmed by summing the intermolecular van der Waals 
energies. That calculation is given in Appendix B. However, zero is not the universal value 
of the Casimir energy for cylinders, as we now remind the reader. 

Of particular interest is the case when £ 2 = 1. With c = 1 in our formulas it corresponds 
to infinitely thin and perfectly conducting cylindrical shell (see Appendix C). Setting £ = 1 
and c = 1 in Eq. ( |3.28| ) we obtain by numerical integration^ 

E = — ^(-0.6517) = -0.05186 4- (4-6) 
47rcr cr 

The sum 2 Y^=\ E n in Eq. ( |3.33| ) can be found by making use of the two leading terms 
in the uniform asymptotic expansion 

^ n 4vra 2 V480 ^ n 2 1792 ^ n 4 J 4vra 2 V480 6 1792 90 J 

= -^0.0210 = 0.0018 4- (4-7) 
47ra 2 a 2 

With higher accuracy (up to 10~ 5 ) this sum has been calculated in [|| by integration of (|3.27|) P| 

00 - 1 1 1 1 

2 E E n - 7T 0-0437 = — — 0.0218 = 0.00174 - . (4.8) 
~[ 4ira z 2 4na z a z 

Substituting Eqs. ( |4.6| ) and ( |4.8| ) into ( ^.33j ) we obtain for the Casimir energy of a perfectly 
conducting cylindrical shell 

E sheU = (-0.1704) = -0.01356 ^ . (4.9) 

This is exactly the result first obtained by DeRaad and Milton |||| . It is worth noting here 
that unlike in that approach the use of the ( function techniques enables us to dispense with 
introduction of the high-frequency cutoff function, although the latter is undoubtedly more 
physical. 

While this paper was being completed the authors became aware of Ref. [27] where the 
vacuum energy of an infinite perfectly conducting cylindrical surface has also been rederived, 
to much higher accuracy, but using a rather more elaborate method. 



2 In the notation of Ref. § this is ~^(S + R + \ ln2vr)— see Eqs. (5.5) and (5.6) of Ref. 
The fn27r term is cancelled by that in Eq. ( |3.33 ) here. 



3 This is exactly the same as — g^-j-R given by Eq. (5.11) of Ref. 
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V. CONCLUSION 



When calculating the Casimir energy there appear to be many arbitrary methods of 
controlling and removing divergences. It is therefore reassuring that unique results emerge 
whatever regularization scheme is adopted. This is somewhat less trivial for this cylindrical 
case than for the case of a sphere, because of the subtleties associated with even space 
dimensions 



28 



However, this universality seems to be characteristic of Casimir calculations 
even in cases where finiteness is not achieved, as for a dielectric ball without the condition 
( p.l| ) imposed [^,^]. Thus our understanding of the Casimir effect seems to be improving. 
However, this comforting conclusion must be tempered by the surprising new result that for 
both dilute dielectric-diamagnetic (satisfying efi = const.) cylinders, and for dilute dielectric 
cylinders, the Casimir energy vanishes. We have as yet no theoretical understanding of these 
zeroes. 
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APPENDIX A: RELATION OF METHOD TO GREEN'S FUNCTION 

APPROACH 

Here we sketch the relation of the formula ( |2.7| ) to the Green's function formalism. We 
will confine our remarks to the case of a massless scalar field, as the generalization to, say, 
electromagnetism is rather immediate. We take the scalar Green's function to satisfy the 
differential equation 

(^-V* ) jG(x,x') = -6(x-x% (Al) 

subject to appropriate boundary conditions. The stress tensor (use of the conformal stress 
tensor has the same effect) is 

1T = d"^d v (j) - i^<9 A 0<9>, (A2) 

so, since the Green's function is given as a vacuum expectation value of a time-ordered 
product of fields, 

G(x,x') = -i{T(j){x)4){x')), (A3) 
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the energy density is 

+00 \ _ 1 r^o qO/ 



u = (t 00 ) = --[d°d ' + V • V']G(x,x') 



(A4) 



To compute the total energy, we integrate u over all space; then we can integrate by parts on 
one of the gradient terms, and use the differential equation (|Al|), omitting the delta function 
because the limit of point coincidence is understood: 

V-V'^-V 2 --^. (A5) 

(The net effect is that the Lagrangian term in t^ u does not contribute.) In terms of the 
Fourier transform of the Green's function 

G(x,x') = [ %L e -**-Vg u {Trf, (A6) 

the Casimir energy is 



(A7) 



Now introduce eigenfunctions of the differential operator V 2 subject to the same bound- 
ary conditions as Q u : 

WW = -k 2 p ^(r), (A8) 

E^(r)^( r ') = 5(r-r'), (A9) 

p 

|(dr)^(r)Vy(r) = V- ( A1 0) 
Then the Green's function has an eigenfunction expansion 

g.(r,r)=E u2 _ kl . (All) 

Carrying out the volume integration, 

/(dr)^(r,r) = £-^, (A12) 

we find that the energy can be written in the form 



OG 



E = -il -e-^-E(^T- + ^T-V (A13) 



Here we have retained a time splitting, r = t — t' — > 0, which is a technique to regulate the 
divergent expression. What does this integral mean? Since the energy must be real, when 
t is set equal to zero, the —i is to be interpreted as an instruction to pick out the negative 
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imaginary part of the integral. That means that the contour of integration C must encircle 
all the poles on the real axis, the positive poles by a contour closed in the counterclockwise 
sense, and the negative poles by a contour closed in the clockwise sense: (see Ref. pTf ) 

E = ~ [ udulnWu- k p )(uj + k p ). (A14) 
ore Jc p 

(We may verify the sign by noting that Eq. ( |2.1|) is formally reproduced.) This is the content 
of 



APPENDIX B: VAN DER WAALS ENERGY OF A DIELECTRIC CYLINDER 

It is now established that for tenuous media the Casimir effect and the sum of molecular 
van der Waals forces are identical [f29j. Here we calculate the latter for a dilute solid cylinder, 
with dielectric constant e 7^ 1 in the interior, e — 1 in the exterior, and fi = 1 everywhere. 



We follow the procedure given in Ref. [ 13]| . The van der Waals energy for this sphere is 

E vdW = -^BN 2 J d D rd D r'{\r ± -v' ± \ 2 + r 2 + r' 2 -2rr'cos9]-~ //2 , (Bl) 

where B = (23/4ir)a 2 , a = (e— l)/47riV being the molecular polarizability, and N being the 
number density of molecules. We have regulated the integral by dimensional continuation, 
D being the number of spatial dimensions, and 7 being the (inverse) power of the Casimir- 
Polder potential. The following calculation is valid providing D > 7; the final result will be 
obtained by violating this condition, by setting D = 3 and 7 = 7. 

We assume translational invariance in the D — 2 transverse directions, so the transverse 
integral is easy: (L is the length of the cylinder and b 2 = r 2 + r' 2 — 2rr' cos 9) 



[ d D ~ 2 r ± d D ~ 2 r' ± [\r ± - r'J 2 + b 2 ]^' 2 = L D ~ 2 f d D ' 2 r ± [r\ + b 2 }^' 2 
J— 00 J — 00 

J-00 Jo t 



r(7/2) 



L 



D-2 



r( 7 /2) Jo 



dtf/ /2-l„-& 



e 



dxe 

00 



-tx* 



D-2 



(L^- 2 (&T /2 - 7/2 - i r(7/2 r + 1} ■ (B2) 



The remaining integral over r, r', 6, 6' is just that given in Ref. [ I3[ . In Eq. (3.25) there, we 
merely set D = 2 and 7 = 7 — D + 2. The result is 

F R/V 2 L °~ 2 2 D -^/ 2 + 1 / 2 r( 7 /2-D/2 + l)r(D/2-7/2 + l/2) 

E vdW - -BN -^- 2 r( 7 /2)r(D/2- 7 /2 + 2)(D- 7 ) ' (B3) 



This is exactly the result found by Romeo []32| . Now when we set D = 3 and 7 = 7 everything 
is finite except for the second gamma function in the denominator, which has a simple pole, 
and thus the Casimir energy vanishes in this case. 
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APPENDIX C: INFINITELY THIN PERFECTLY CONDUCTING 

CYLINDRICAL SHELL 



We show here that the Casimir energy for an infinitely thin perfectly conducting cylin- 
drical shell is given by Eqs. ( |3.18|) and (|3.19| ) with £ 2 = 1. In this for a perfectly 
conducting spherical shell, the frequencies of electromagnetic oscillations inside and outside 
the shell turn out to be different [16]. The frequencies of the TE- modes are determined byQ 



For the TM-modes we have 



In these equations 



J' n (\a) 
H' n (Xa) 



J n (Xa) 
H n (\a) 



0, r < a, 
0, r > a . 



0, r < a, 
0, r > a . 



n = 0,±1,±2, 



(CI) 
(C2) 



(C3) 
(C4) 



(C5) 



Substituting A™ and A™ in Eq. ( |3.9D by the new equations (PT|)-(|C4]) we obtain [cf. with 
Eq- ( TO ? 



E 



Ana 2 



E 



>dy In {-Ay 2 I n (y)I' n (y)K n (y)K' n 



(C6) 



To rearrange the argument of the logarithmic function in this formula we again use equalities 
( gig ) and fl3l7D . This gives Eqs. ( gjjp and (|3T9|) with ^ 2 = 1 and c = 1: 

E = ^J2 I ydy\n{l-[y(I n (y)K n (y))f} , (C7) 

which is, of course, the unregulated version of the result derived rigorously in Ref. B. 



4 The exterior modes can be considered only formally, since the Hankel functions have only a finite 
number of complex zeroes. See Ref. Q. Nevertheless, this formal procedure yields the correct 
result. The mode sum breaks down in this case, because the singularity structure is not that 
assumed in ( p.7| ), but the Green's function method retains meaning. 
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